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ABSTRACT 

In this work we analyse late-time (t > 100 d) optical spectra of low-redshift (z < 0.1) 
Type la supernovae (SNe la) which come mostly from the Berkeley Supernova la 
Program dataset. We also present spectra of SN 201 lby for the first time. The BSNIP 
sample studied consists of 34 SNe la with 60 nebular spectra, to which we add nebular 
spectral feature measurements of 20 SNe la from previously published work (Maeda 
et al. 2011; Blondin et al. 2012), representing the largest set of late-time SN la spectra 
ever analysed. The full width at half-maximum intensity (FWHM) and velocities of the 
[Fe III] A4701, [Fe II] A7155, and [Ni II] A7378 emission features are measured in most 
observations of spectroscopically normal objects where the data have signal-to-noise 
ratios > 20 px _1 and are older than 160 d past maximum brightness. The velocities 
of all three features arc seen to be relatively constant with time, increasing only a few 
to ~20 km s _1 d _1 . The nebular velocity (u ne b, calculated by taking the average of 
the [Fe II] A7155 and [Ni II] A7378 velocities) is correlated with the near-maximum- 
brightness velocity gradient and early-time ejecta velocity. Nearly all high velocity 
gradient objects have redshifted nebular lines while most low velocity gradient objects 
have blueshifted nebular lines. No correlation is found between v nc h and Ami5(£>), 
and for a given light-curve shape there is a large range of observed nebular velocities. 
The data also indicate a correlation between observed (B — V) and w n ob- 
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1 INTRODUCTION 

Type la supernovae (SNe la) can be used to accurately mea- 



sure cosmological parameters (e.g., Astier et al.|2006 


Riess 


et al. 2007; 


Wood-Vasey et al.||2007| |Hicken et al. 


2009 


Kessler et a 


. 2009 


Amanullah et al.||2010| |Conley et al.| 


2011 Suzuki et al. 


20121, and they led to the discovery of 



the accelerating expansion of the Universe ( Riess et al.|1998 



Perlmutter et al.|1999[ ) . SNe la are thought to be the result of 
thermonuclear explosions of C/O white dwarfs (WDs) (e.g., 
Hoyle fc Fowler|1960||Colgate fc McKee|1969||Nomoto et aT] 



At these late epochs, SN la spectra consist of broad emis- 
sion lines of (mostly) iron-group elements (IGEs) and can 
yield valuable insights into the physics of the explosion it- 
self. However, SNe la often appear quite faint at these late 
phases and thus not many late-time spectra exist in the lit- 
erature. 

There are a handful of relatively "normal" SNe la (i.e., 
ones that follow the light-curve width versus luminosity re- 
lationship; Phillips 1993 1 that have published spectra at 
phases later than ~ 100 d (e.g., Stritzinger & Sollerman| 



1984| see |Hillebrandt fc Niemeyer|2000| for a review). Despite |2007| |Stanishev et al.||2007| |Leloudas et al.||2009| >. In ad- 



their cosmological utility, we are still missing a detailed un- 
derstanding of the progenitor systems and explosion mech- 
anisms (sec Howell 2011, for further information). 

By about 100 d past maximum brightness, SN la ejecta 
have expanded significantly and the SN enters the so-called 
"nebular phase" (as opposed to the "early-time" phase). 



E-mail: jsilverman@astro.as.utexas.edu 



dition, there are some rare and peculiar SNe la with pub- 
lished late-time spectra (e.g., |Jha et al.|20 06 Silve rman et al. 
2011 1. While there are large spectroscopic samples of SNe la 
published jMatheson et al.|2008||Blondin et al.|2012| [Silver- 



man et al.|2012a[ ), only a tiny fraction of those data consist 
of late-time spectra. However, moderate-sized comparative 
studies of ~ 14-24 SNe la with late-time spectra have been 
undertaken ( jMazzali et al.|1998|[Maeda et al.|2010a||Blondin 
et al.|2012| |. 
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Lick Observatory with the Kast double spectrograph ( Miller 



4000 5000 6000 7000 8000 9000 10000 
Rest Wavelength (A) 

Figure 1. Spectrum of SN 2011by from 206 d past maximum 
brightness with prominent features labeled. The left-hand shaded 
area represents the [Fe III] A4701 feature while the right-hand 
shaded area represents the [Fe II] A7155 and [Ni II] A7378 fea- 
tures. 



These works have concentrated mainly on three broad 
emission features centred near 4701, 7155, and 7378 A ( see 
Figure [IJ. The 4701 A feature is likely a blend of various 
[Fe III] lines ( |Stritzinger et alpOOo] |Maeda et al.||2010b l 
which come from material produced by a supersonic burn- 
ing front (i.e., a detonation; Maeda et al. 2010a I. On the 
other hand, the 7155 A feature is likely due to [Fe II] ( Maeda 



et al.|2010a| and the 7378 A feature is probably from [Ni II 
(Maeda et al. 2010b). These are both thought to trace mate- 
rial from where the explosion began and where the burning 
initially proceeded subsonically (i.e., a deflagration; |Maeda| 
et al.|2010a l. 

To the previous late-time SN la spectral studies we 
now add this work, where we analyse 60 nebular-phase 
(t > 100 d) low-resolution optical spectra of 34 low-redshift 
(z < 0.1) SNe la obtained as part of the Berkeley SN la 
Program (B SNIP; |Silverman et al.||2012a] ). The data are 
presented by Silverman et al. (2012a I (supplemented herein 
with one new object: SN 2011by)[] and we utilise spectral 
feature measurement tools similar to those described by|Sil-| 



verman et al. (2012dl. The spectral data used in this work 



are summarised in Section[2] and our method of spectral fea- 
ture measurement and our determination of nebular veloci- 
ties are described in Section [3] Comparison of our late-time 
spectral measurements to other SN la observables can be 
found in Section [4] We present our conclusions in Section [5] 



& Stone||1993 1 and the 10 m Keck telescopes with the Low 



Resolution Imaging Spectrometer (LRIS; Oke et al. 19951. 
The Kast data typically cover 3300-10,400 A with resolu- 
tions of ~11 and ~6 A on the red and blue sides (crossover 
wavelength ~5500 A), respectively. The LRIS spectra usu- 
ally have a range of 3100-9200 A with resolutions of ~7 
and ~5.5 A on the red and blue sides (crossover wavelength 
~5600 A), respectively. 

All data were reduced using standard reduction meth- 
ods. For more information regarding the observations and 



data reduction, see [Silverman et al. (2012a I. The spectral 
ages of the BSNIP data referred to throughout this work are 
calculated using the redshift and Julian Date of maximum 
presented in Table 1 of Silverman et al. (2012a I. Further- 



more, photometric parameters (such as light-curve width 
and colour information) used in the present study can be 
found in Ganeshalingam et al. (20101. 



The BSNIP sample contains 81 spectra of 43 SNe la 
older than 100 d past maximum brightness. In order to 
measure the three emission features mentioned above (i.e., 
[Fe III] A4701, [Fe II] A7155, and [Ni II] A7378), we required 
that the spectra span ~4300-7250 A, which decreased the 
sample to 73 spectra of 38 SNe la. In previous BSNIP studies 
( Silverman et al.|2012d I, we a priori ignored the extremely 
peculiar SN 2000cx (e.g.,|Li et al.|2001[ ), SN 2002cx (e.g.,p 



et al |2003||Jha et al.|2006| ), SN 2005gj (e.g., |Aldering et al 



2006 Prieto et al. 20071, and SN 2005hk (e.g., Chornock 



et al. 2006 Phillips et al. 2007 1. These objects are so spectro 



scopically peculiar that it is difficult to consistently measure 
their spectral features in comparison to the bulk of the SN la 
population. Thus, in this work we will only concentrate on 
the SNe la that follow the "Phillips relation" ( |Phillips|1993 \ 
and can be used as cosmological distance indicators. We also 
include the underluminous SN 1991bg-like objects, which lie 
on a quadratic relationship between light-curve width and 
luminosity (e.g., Filippenko et al.||1992b Leibundgut et al 



1993) 



After removing the peculiar objects mentioned above, 
we are left with 33 SNe la with 58 BSNIP spectra having 
t > 100 d. Of these, 7 were in the late-time sample studied 
by |Maeda et al.|(|2010a[), and those 7 plus 2 more were in the 
sample of Blondin et al. (2012). Each of the objects studied 



herein and the spectral phases are listed in Table [T] along 
with the full width at half-maximum intensity (FWHM) and 
velocity of the [Fe III] A4701 feature, and the velocities of 
the [Fe II] A7155 and [Ni II] A7378 features. 



2 DATASET 



The late-time spectral sample used in this work is a subset 
of spectra presented by Silverman et al. (2012a I, with the 
addition of one new object (see Section |2.1[ ). The majority of 
the spectra were obtained using the Shane 3 m telescope at 



1 All of the BSNIP data, along with the SN 201 lby data, will be 
made public on 1 Jan. 2013 via the SuperNova DataBase (SNDB; 
Silverman et al.|2012ai. 
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2.1 SN 2011by 

To the BSNIP sample discussed above we also add new data 
for the nearby Type la SN 2011by (see Appendix [A|for more 
information). We obtained broad-band BVRI photometry 
as well as 7 near-maximum-brightness spectra and 2 late- 
time spectra of SN 2011by. All of these data indicate that 
SN 2011by is photometrically (Am 15 (B) = 1.14±0.03 mag) 
and spectroscopically normal. 

SN 2011by will serve as a high-quality individual case 
study that can be directly compared to the larger late-time 
BSNIP sample discussed above. Adding this object and its 
two late-time spectra to the aforementioned dataset yields 
a sample of 60 spectra of 34 SNe la with t > 100 d past 
maximum brightness. This represents one of the largest sets 
of late-time SN la spectra ever analysed. 



3 MEASURING NEBULAR SPECTRAL 
FEATURES 

The algorithm used in this work to measure the [Fe III] 
A4701, [Fe II] A7155, and [Ni II] A7378 emission features 
is similar to that used to measure absorption features in 
near-maximum-brightness BSNIP spectra. It is described in 



detail by Silverman et al. ( 2012d[ ), but here we give a brief 
summary of the procedure. Each spectrum first has its host- 
galaxy recession velocity removed and is corrected for Galac- 
tic reddening (according to the values presented in Table 1 
of Silverman et al. 2012a I, and then is smoothed using a 



Savitzky-Golay smoothing filter ( Savitzky fc Golay|1964 l. 



For each feature, a linear background is denned by con- 
necting points on either side of the feature and then sub- 
tracting out that background flux level. A cubic spline is 
then fit to the smoothed data and the expansion velocity 
is calculated from the wavelength at which the spline fit 
reaches its maximum. This maximum is also used to mea- 
sure the FWHM for the [Fe III] A4701 feature. 

[Fe III] A4701 is usually the strongest feature in late- 
time SN la spectra and is often seen to be a single-peaked, 
distinct feature. On the other hand, [Fe II] A7155 and [Ni II] 
A7378 are found to be blended into either a single- or double- 
peaked broad emission feature. Any spectrum where the two 
features were severely blended with each other (or showed 
more than two peaks near the wavelengths of interest) did 
not have velocities measured for [Fe II] A7155 and [Ni II] 
A7378. For the spectra that do show clearly double-peaked 
profiles, a cubic spline is again fit to the smoothed data and 
the expansion velocities are calculated from the wavelengths 
at which the spline fit reaches its maximum on either side of 
the local minimum (i.e., the point between the two peaks). 
The FWHM and velocity of the [Fe III] A4701 feature and 
the velocities of the [Fe II] A7155 and [Ni II] A7378 features 
are displayed in Table [l] 

While we were able to measure [Fe III] A4701 in all 
late-time spectra presented herein whose wavelength range 
included this feature, we were only able to accurately mea- 
sure the velocities of [Fe II] A7155 and [Ni II] A7378 in 22 
spectra of 15 SNe la (nebular velocities of 6 of these ob- 
jects have been calculated previously; Mae da et al.||2010a[ 
|Blondin et al.|[2012[ ). In these previous studies, velocities 
were reported even when there was only a single peak near 



these wavelengths. With only one peak, however, it is not 
clear which reference wavelength to use when calculating a 
nebular velocity, and so for a velocity measurement in this 
work we require that exactly two peaks be present in a late- 
time spectrum. 



The underluminous SN 1991bg-like objects (e.g., Filip- 
|penko et al.||1992b Leibundgut et aT7||1993 1 have narrower 



emission profiles than do normal SNe la at late times, but 
they do not typically show the [Fe II] A7155 and [Ni II] A7378 



features. Furthermore, the SN 1991T-like objects (e.g., Filip 
penko et al.||1992a| [Phillips et al.|1992|> and SN 1999aa-like 



objects (e.g., Li et al. 2001 Strolger et al.|[2002 Garavini 



et al. 20041, which together generally represent overlumi- 



nous SNe la, often have complex and multi-peaked profiles 
near 7300 A, and thus we are again unable to accurately 
measure the velocities of [Fe II] A7155 and [Ni II] A7378 
for most of these objects. We are, however, able to measure 
velocities in one spectrum of SN 1999aa. Note that Maeda 
et al. (2010a I cut these peculiar objects from their sample 



as well. Of the remaining "normal" SNe la, 6 of the spectra 
in which we are unable to reliably measure these features 
have low signal-to- noise ratios (S/N < 20 px _1 ), with 3 of 
those having S/N < 5 px _1 , and all normal SN la spectra 
with S/N > 20 px _1 were successfully measured. 

Even after removing 9 spectra of peculiar SNe la and 
6 spectra with low S/N, there are still 23 instances where 
we are unable to measure reliable velocities for the [Fe II] 
A7155 and [Ni II] A7378 features. Nearly all of these failures 
are due to the spectra being obtained at too young a phase. 
Figure [2] shows a histogram of the ages of the spectra in- 
vestigated in this work, after removing peculiar objects and 
low-S/N spectra. The grey histogram represents the ages of 
spectra in which we are unable to measure [Fe II] A7155 and 
[Ni II] A7378 velocities, while the white histogram represents 
ages where we successfully measure the velocities. At epochs 
less than ~ 160 d past maximum brightness, the emission 
near 7300 A is often still multi-peaked and too complex to 
definitively identify the [Fe II] A7155 and [Ni II] A7378 fea- 
tures. Thus, in order to be able to measure these velocities, 
one should wait until ~ 160 d past maximum, after which 
time we are able to measure the velocity in 16 of 17 spectra. 
Despite this finding, there are some BSNIP spectra at 115, 
118, 126, 147, and 149 d past maximum that have measured 
velocities. 

The oldest spectrum in our sample is from 362 d past 
maximum brightness, but it is SN 1991bg-like and we are un- 
able to measure the velocities of the [Fe II] A7155 and [Ni II] 
A7378 features. The next-oldest spectrum is from 358 d past 
maximum and in that case we are able to measure the veloc- 
ities accurately. Thus, the epoch range in which we success- 
fully measure nebular velocities of [Fe II] A7155 and [Ni II] 
A7378 is 115-358 d past maximum, with the majority of the 
spectra in the range 160-358 d past maximum. Both |Maeda| 
|et al.| ( |2010a[ ) and |Blondin et al.| ( |2012[ ) had a similar range 
(^150-400 d past maximum), with the vast majority of their 
spectra being older than ~200 d past maximum. 

3.1 Calculating Nebular Velocities 



Following Maeda et al. (2010a I and Blondin et al. (20121, 



we calculate nebular velocities from the mean of the individ- 
ual velocities of the [Fe II] A7155 and [Ni II] A7378 features. 
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is approximately ±2900 km s , similar to the range found 
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Figure 2. Histogram of spectral ages after removing peculiar ob- 
jects and low-S/N spectra. The grey area represents spectra in 
which we are unable to measure [Fe II] A7155 and [Ni II] A7378 
velocities. The white area represents spectra in which we can mea- 
sure these velocities. Note that the vast majority of spectra with 
t > 160 d have measured velocities. 



Thus, the nebular velocity represents the typical velocity of 
material synthesized at the earliest times in the SN explo- 
sion. Previous studies used the difference between the two 
velocities as the uncertainty of the nebular velocity, but we 
use the difference divided by two (since that is how much 
the average differs from each measurement). 

For the 4 SNe la in which we are able to measure both 
velocities in multiple late-time spectra, the nebular velocity 
is the average of all [Fe II] A7155 and [Ni II] A7378 velocities 
and the uncertainty is the average difference. The calculated 
nebular velocities (v ae b) of all 15 SNe la with [Fe II] A7155 
and [Ni II] A7378 velocity measurements are listed in the 
ninth column of Table [2] along with various classifications 
and other measured properties. 

The second column of Table [2] gives the decline rate of 
the light curve, Amis(i?), and the third column shows the 
observed B — V colour of the SN at B-band maximum bright- 
ness, (B — V) , both taken from IGaneshalingam et al 



( 2010 1. Also displayed in the table is the velocity of the Si II 



A6355 feature within 8 d of maximum brightness (fifth col- 
umn) and a classification based on that velocity of either 
normal velocity or high velocity (fourth column; e.g., Wang 



et al.|2009 I . The rate of decrease of the expansion velocity, 
also known as the velocity gradient (t\ seventh column), and 
the velocity of the Si II A6355 feature extrapolated to t = d 
using v (vo, eighth column) are also shown in Table |2| along 
with a classification based on v (sixth column; e.g., |Benetti 
et al. |2005" \ . The tenth (eleventh) column of the table is the 
EW of narrow Na I D (Ca II H&K) absorption from the 
host galaxy of each SN (see Section 4.3.1 for more informa- 
tion on these measurements). Objects listed with a 'f are 
part of the low-extinction sample (see Section 4.3.21. Note 



that while most of the measurements come directly from 
Silverman et al. ( 2012d I , some of the near-maximum-light 



Si II A6355 velocities were measured from data presented by 
Matheson et al.l (12008b andlBlondin et all (12012 I. 



in previous work (Maeda et al. 2010a Blondin et al. 2012) 



A comparison of the nebular velocities of the 6 SNe la that 
are found in the current sample as well as previous stud- 
ies shows that the values calculated herein are consistent 
with previous work, with average absolute differences of a 



few hundred km s and no detectable bias (Maeda et al 



2010a||2011||Blondin et al.||2012| >. The average uncertainty 



of Vneb (which is half the difference between the individual 
velocities of the [Fe II] A7155 and [Ni II] A7378 features) is 
~490 km s" 1 , slightly larger than the typical uncertainty in 
the measurements of Maeda e t al.| (2010a I. 



4 ANALYSIS 



4.1 [Fe III] A4701 



The FWHM of the [Fe III] A4701 feature has been seen 



to possibly be anticorrelated with Amis(_B) (Mazzali et al. 
|1998[ ). In that work, both overluminous SN 1991T-like ob- 
jects and underluminous SN 1991bg-like objects were found 



to follow the correlation. More recently, Blondin et al. ( 2012 1 



found a similar result with a Pearson correlation coefficient 
of —0.71 for all of their late-time data, but when their two 
underluminous SNe la were removed it decreased to —0.17 
(i.e., no correlation). This latter result was questioned, how- 
ever, due to the fact that spectra at epochs earlier than 
250-300 d past maximum brightness were used (Mazzali & 
Hachinger||2012 l. 



Figure |3| shows the FWHM measurements of the [Fe III] 
A4701 feature versus Ami 5 (B) from the BSNIP spectra, as 
well as previously published values for 18 SNe ( jMaeda et al.| 



2011 Blondin et al. 20121. The shape of each data point 
signifies its "Benetti Type" (based on v): squares are low 
velocity gradient (LVG) objects, triangles are high veloc- 
ity gradient (HVG) objects, stars are FAINT objects, and 
circles are objects for which we are unable to determine a 
Benetti Type. The colour of each point represents its "Wang 
Type" (based on Si II A6355 velocity): blue is for normal- 
velocity objects, red is for HV objects, and black is objects 
for which we are unable to determine a Wang Type. Some 
of the outlying points are labeled with their object name. 



As in previous studies (Mazzali et al. 1998| |Blondin 



|et al.|2012| ), the BSNIP data show that the FWHM of [Fe III 
A4701 is anticorrelated with Ami5(B) (Pearson correlation 
coefficient of —0.66) and the linear fit to all of the data 
is shown as the solid line in Figure [3] The slope of this 
line, —0.073 mag/(10 3 km s _1 ), is consistent with that of 



Blondin et al. (2012) but inconsistent with that of Mazzali 



et al.| ( |1998[ ), and the correlation is significant at the > So- 



level. If we remove only the underluminous SN 1991bg, the 
linear fit is nearly unchanged (the dashed line in the figure) 
and the correlation weakens only slightly (Pearson coeffi- 
cient of —0.61). If we instead remove the 8 objects with 
Ami5(_B) > 1.6 mag (i.e., all of the underluminous objects 
in the dataset), the correlation effectively disappears (Pear- 
son coefficient of —0.34) and the slope of the linear fit (the 
dotted line in Figure [3|) is consistent with at the ~2<r level. 



3 |3| 

r[Bi 



The findings of Blondin et al. (20121 confirmed above 



The range of nebular velocities calculated in this study 



have recently been questioned due to their use of spectra 
at epochs earlier than 250-300 d past maximum brightness 
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Table 2. Summary of objects studied. Uncertainties are given in parentheses. 



SN Name 


Ami 5 (S) a 


(S - v) max a 


Wang 


Si II A6355 Vel.° 


Bcnctti 






«0* 


Nebular Vcl. g 


EW(Na) 


EW(Ca) 




(mag) 


(mag) 


Typc b 


(10 3 km s _1 ) 


Typc d 


(km s 


-' cl- 1 ) 


(10 3 km s" 1 ) 


(km s 




(A) 


(A) 


SN 1989M 


1.10 (0.20) 




HV 


12.46 (0.10) 


HVG 


291.84 


(139.49) 


13.19 (0.42) 


560 


(110) 


1.88 (0.73) 


0.36 (0.18) 


SN 1993Z 


















1990 


(630) 


1.56 (0.40) 


0.96 (1.67) 


SN 1994D 1 " 


1.36 (0.05) 


-0.053 (0.035) 


N 


10.72 (0.10) 


LVG 


33.43 


(6.77) 


10.60 (0.05) 


-1720 


(570) 


0.40 (0.44) 


0.09 (0.08) 


SN 1998bu 


1.05 (0.03) 


0.292 (0.032) 


N 


10.84 (0.11) h 


LVG h 


46.52 


(5.38) h 


11.01 (0.05) h 


-1800 


(280) 


1.52 (1.12) 


0.39 (0.19) 


SN 1999aat 


0.79 (0.05) 


-0.093 (0.033) 




10.50 (0.10) 


LVG' 


14.33 


(3.52)' 


10.49 (0.03)' 


760 


(70) 


0.44 (0.64) 


0.27 (0.14) 


SN 1999ght 


1.69 (0.05) 


0.190 (0.000) 


N 


11.02 (0.10) 


FAINT 


46.49 


(10.13) 


11.20 (0.13) 


1570 


(280) 


0.50 (0.25) 




SN 2002bo 


1.15 (0.04) 


0.315 (0.052) 


HV 


13.88 (0.10) 


HVG 


245.14 


(8.12) 


13.61 (0.09) 


2430 (100) 


2.36 (0.53) 


0.57 (0.31) 


SN 2002fkt 


1.20 (0.03) 


-0.155 (0.032) 


N 


9.51 (0.10) 


LVG' 


11.21 


(8.28)' 


9.74 (0.05)' 


-2060 


(690) 


0.46 (0.09) 


0.26 (0.07) 


SN 2003gs 


1.83 (0.02) k 


0.634 (0.021) k 




11.00 (0.20) k 










2540 (1330) 


0.98 (0.49) 




SN 2004dt 


1.29 (0.05) 


-0.080 (0.034) 


HV 


14.34 (0.10) 


HVG 


269.45 


(8.34) 


14.70 (0.10) 


-2850 (1380) 


1.80 (1.04) 




SN 2005cft 


1.08 (0.03) 


-0.062 (0.031) 


N 


10.26 (0.10) 


HVG 


106.69 


(149.59) 


10.13 (0.25) 


910 


(500) 


0.14 (0.09) 


0.53 (0.26) 


SN 2006X 


1.10 (0.04) 


1.260 (0.034) 


HV 


15.26 (0.10) 


HVG' 


184.80 


(5.82)' 


15.63 (0.06)' 


2700 


(410) 


2.04 (1.61) 


1.05 (0.46) 


SN 20071e 


1.02 (0.04) 


0.280 (0.039) 


HV 


12.09 (0.10) 


HVG 


93.35 


(12.65) 


12.77 (0.18) 


2030 


(140) 


2.36 (0.62) 


0.33 (0.17) 


SN 2008Qt 


1.25 (0.08) 


-0.012 (0.031) 


HV 


11.09 (0.10) 


HVG' 


82.66 


(16.22)' 


11.64 (0.05)' 


-2200 


(750) 


0.44 (0.22) 




SN 2011by 


1.14 (0.03) 


-0.061 (0.032) 


N 


10.35 (0.14) 


LVG 


52.83 


(8.76) 


10.40 (0.08) 


-1450 


(170) 


0.66 (0.16) 


0.47 (0.41) 



high velocity; 'N' — normal. 



^Object is part of t he low- extinct ion sample (se e text for details). 
a Values taken from |Ganeshalingam et al.j|2010j l. 

k Classification based on the velocity of the Si II A6355 line near maximum brightness jWang et al. 20091. 'HV 
Taken from [Silverman et al.|l2012dl. 

c Velocity or the Si 11 At>355 feature within 8 d of maximum brightnes s. Taken from [Silverman et al.j l j2012d| unless otherwise noted. 
^Classification based on the velocity gradient of the S i II A6355 line ^Benetti et al.|2005| . 'HVG' high velocity gradient; 'LVG' — low velocity 
gradient; 'FAINT' — faint/undcrluminous. Taken from |Sil verm an et al.||2012d^ unless otherwise noted. 

B Change in velocity with time, near maximum brightness (i.e., the velocity gradient). Taken from |Sjlyerman et alH j2012d| unless otherwise noted. 
^Velocity of the Si II A6355 feature extrapolated to t — d based on the velocity gradient. Taken from Silverman ct al. ( 2012d I unless otherwise 
noted. 

g Average nebular velocity of the [Fe II] A7155 and [Ni II] A7378 features for all latc-timc spectra of a given object. The uncertainty is the average 

difference between the velocity of the two features. 

h Si II A6355 velocity. Bcnctti type, velocity gradient, and vq measured from data presented by |lVIatheson et al.| l j2008| . 
Classification scheme only applies to spcctroscopically normal SNe la. 

J Bcnctti type, velocity gradient, and vq measured from BSNIP d ata as well as data presented by 
' ! Ami5(B) and (B — V r ) max taken fromlKrisciunas et al.l a2009l . Si II A6355 velocity taken from 



Blondin et al. 


2012[ 


Suntzeff et al. 


2003 1 



2.0 
1.8 

m 1.6 

S. 

g 1.4 
I 1.2 
1.0 
0.8 




(Mazzali & Hachinger 20121. To investigate this possibil- 
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FWHM (10 3 km s" 1 ) of [Fe III] >.4701 



Figure 3. FWHM of [Fe III] A4701 ver sus Am 15 (B) from BSNIP 
and 18 objects from |Maeda et al.lpoTi] and |Blondin et a.1. | ( [2012] ) . 
Squares are LVG objects, triangles are HVG objects, stars are 
FAINT objects, and circles are objects for which we are unable 
to determine a Benetti Type. Blue points are normal- velocity ob- 
jects, red points are HV objects, and black points are objects for 
which we are unable to determine a Wang Type. Some of the 
outlying points are labeled. Linear fits to the data are also dis- 
played: all data (solid), after removing SN 1991bg (dashed), and 
after removing all 6 objects with Amirj(B) > 1.6 mag (dotted). 



ity, we performed a separate analysis of the FWHM of the 
[Fe III] A4701 feature versus Amis(B) using only the oldest 
spectra from BSNIP, Maeda et al. (2011 1, and Blondin et al. 



(2012|. There were 16 (8) objects with spectra obtained at 
t > 250 d (300 d) past maximum brightness, and only a mod- 
est anticorrelation was found for these data with a Pearson 
coefficient of —0.55 (—0.69). When removing the 2 (1) under- 
luminous objects with Amis(B) > 1.6 mag having spectra 
at these latest epochs, the possible correlation completely 
disappears. Thus, it seems that at all epochs the possible re- 
lationship between the FWHM of the [Fe III] A4701 feature 
and Ami5(B) is being driven almost completely by the most 
underluminous SNe la, perhaps indicating the existence of 
two groups of points rather than a linear correlation. 

Regarding the velocity of the [Fe III] A4701 feature, 



Stritzinger et al. ( 2006 1 found no significant temporal change 



in their measurements and Maeda et al 
the feature shows 



( 2010a[ ) state that 
However, 



'virtually no Doppler shift.' 
later work from the latter group show blueshifts of 2000- 
4500 km s" 1 for t < 200 d, with the velocity steadily ap- 
proaching for t > 200 d JMaeda et al.||2010bl 120111). This 



more recent result is also seen in the BSNIP data. In Fig- 
ure [4] all measured velocities of [Fe III] A4701 are plotted, 
with solid lines connecting multiple velocities of individual 
objects. Squares represent spectra for which we measured 
a nebular velocity, while x symbols represent spectra for 
which we did not. 

The overall trend (which matches what was found by 
Maeda et al.||2010b l is clear: velocities are blueshifted typi- 
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Figure 4. All velocities of [Fe III] A4701 measured in this 
work. Solid lines connect multiple velocities of individual objects. 
Squares represent spectra for which we measured a nebular ve- 
locity; X symbols represent spectra for which we did not. 



cally by ~2000-4000 km s _1 at t = 100 d and the amount of 
blueshift decreases approximately linearly with time. This 
temporal evolution is quite slow, however, with a typical 
change in velocity of only ~10-20 km s _1 d _1 . The velocity 
of the [Fe III] A4701 feature is uncorrelated with the nebu- 
lar velocity and hence with the individual velocity of either 
[Fe II] A7155 or [Ni II] A7378. Furthermore, the [Fe III] A4701 
velocity is also uncorrelated with the Si II A6355 velocity and 



4.2 [Fe II] A7155 and [Ni II] A7378 

The individual velocities of the [Fe II] A7155 and [Ni II] 
A7378 features show a relatively large scatter at all epochs, 
which is reflected in the large range of calculated nebular 
velocities. The temporal velocity evolution of both features 
is even slower than that of the [Fe III] A4701 feature, with 
typical changes of a few km s _1 d _1 . 

Since the average of the velocities of [Fe II] A7155 and 
[Ni II] A7378 is used to calculate fj nc b, it is instructive to 
compare these two velocities with each other. The differ- 
ence in these velocities has a relatively large scatter at all 
epochs, but we find no bias above or below zero velocity. 
The average absolute difference between the two velocities, 
however, is ~900 km s , which is reflected in our v nc ^, un- 
certainties discussed above. Despite this, the two velocities 
are highly correlated (Pearson coefficient of 0.84) and the 
difference between the two has a typical temporal evolution 
of < 10 km s _1 d _1 . 



4.3 Nebular Velocities 

4-3.1 Spectral Comparisons 

Like the [Fe II] A7155 and [Ni II] A7378 velocities from which 
it is calculated, the nebular velocity does not change much 
with time (typical values being only a few km s _1 d _1 , sim- 
ilar to what has been seen previously; e.g., |Maeda et al.| 
20111. Furthermore, all four SNe la which have measured 



Figure 5. Nebular velocity versus velocity gradient. Shapes and 
colours of the data points are the same as in Figure [3] The ver- 
tical dashed line at km s _1 separates blueshifted (left) and 
redshifted (right) nebular velocities. The horizontal dashed line 
at 70 km s — 1 d — 1 separates LVG (below) and HVG (above) ob- 
jects. 



[Fe II] A7155 and [Ni II] A7378 velocities in multiple late- 
time spectra show no significant change in their nebular ve- 
locities with time (for 130 < t < 360 d). Thus, a single 
measurement of v nB b appears to be sufficient to determine 
the typical nebular velocity of a given SN la. 

Measurements of the velocity gradient were made for 
all objects with multiple near-maximum-brightness spectra 
in BS NIP, in |Matheson et al 



Blondin et al 

( 2012t > (i-e., all SNe la listed in Table|2|except SNe 1993Z and 



2003gs). The v values calculated in this work are consistent 



at the 2o level with those found in previous studies ( Benetti 



et al. 2005 Maeda et al. 2010a I . We also plot previously 



published v and v ne h values for 12 SNe la from |Maeda et aL] 
( |2011| ) and |Blondin et al.| ( [2012 1. 

Of the 25 objects with measured velocity gradients, 11 
are LVG, 12 are HVG, and 2 are FAINT, which is sim- 
ilar to the distribution of types found by |Maeda et al.| 
| |2010a| ), though they have a slightly larger fraction of LVG 
objects. Figure [5] presents these velocity-gradient measure- 
ments versus the nebular velocity; shapes and colours of 
the data points represent the Benetti and Wang types, 
respectively, as in Figure [3] The vertical dashed line at 
km s _1 is the demarcation between blueshifted (left) and 
redshifted (right) nebular velocities; the horizontal dashed 
line at 70 km s _1 d _1 is the demarcation between LVG (be- 
low) and HVG (above) objects. 

Maeda et al.| ( |2010a| l found a correlation between n nc b 



and i), and the BSNIP data support this conclusion, with 
some caveats. The Pearson coefficient for all of the data in 
Figure [5] is only 0.25 and the Spearman rank correlation 
coefficient is 0.390 However, SN 2004dt was found to be 
somewhat peculiar (showing evidence for high-velocity and 
clumpy ejecta at early times; e.g., |Altavilla et al.|2007| ) and 



2 This is perhaps the more appropriate value to examine here 
since the relationship appears to be nonlinear. 
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Nebular Velocity (10 3 km s" 1 ) 



Figure 6. Nebular velocity versus Si II A6355 velocity. Shapes 
and colours of the data points are the same as in Figure [3] Some 
of the outlying points are labeled. The vertical dashed line at 
km s _1 separates blueshifted (left) and redshiftcd (right) neb- 
ular velocities. 



did not follow the relationship of Maeda et al. ( 2010a I. Also, 



there are only two BSNIP spectra of SN 1989M and they 
were obtained just 1 d apart, thus the lever arm for the v 
calculation is small and the accuracy of the measurement is 
small (hence the large error bars for that data point). 

Removing SN 2004dt (SNe 2004dt and 1989M) yields 
a Pearson correlation coefficient of 0.51 (0.62) and a Spear- 
man rank correlation coefficient of 0.56 (0.59). This indicates 
that the nebular velocity is correlated with the velocity gra- 
dient, with the significance of the correlation at the 2a (3a) 
level. If we ignore SN 2004dt, all HVG objects have red- 
shifted nebular lines, except SN 2008Q which is extremely 
close to the LVG/HVG border. We also find that nearly 3/4 
of all LVG objects have blueshifted nebular spectral features. 
A Kolmogorov-Smirnov (KS) test on the nebular velocities 
of LVG and HVG objects implies that they very likely come 
from different parent populations (p ss 0.016). When remov- 
ing SN 2004dt, the difference is even stronger (p « 0.007). 

The near-maximum-brightness Si II A6355 velocity can 
also be compared to the nebular velocity for 14 of the SNe la 
in Table [2} as well as 15 SNe la from |Maeda et al.| ( pf3TT] > and 



Blondin et al. (2012 1. We measure the Si II A6355 velocit 



only in BSNIP spectra within 8 d of maximum brightness 
When multiple BSNIP spectra within 8 d of maximum exist 
for a given object, only the one closest to maximum bright- 
ness is used for the Si II A6355 velocity measurement. The 
near-maximum-light velocities are plotted against the neb- 
ular velocities in Figure [6] with shapes and colours of data 
points the same as in Figure |3j the vertical dashed line sepa- 
rates blueshifted (left) and redshifted (right) nebular veloc- 
ities. 

A moderate correlation exists for all of the data pre- 



3 Spectra within 5 d of maximum brightness were used by |Sil-| 
jverman et al.| | |2012d[ l , but applying this constraint would remove 
three objects from the study of Si II A6355 velocity versus nebular 
velocity. 



sented in Figure [fj] (Pearson and Spearman coefficients are 
both 0.51). If SN 2004dt is again ignored, a strong correla- 
tion exists (Pearson coefficient of 0.70 and Spearman coef- 
ficient of 0.64). Most HV objects have redshifted emission 
lines at late time, but normal-velocity objects have both 
blueshifted and redshifted nebular features. A KS test of 
the nebular velocities indicates, however, that the HV and 
normal-velocity objects do come from different parent pop- 
ulations (p ~ 0.021 for all of the data and p ~ 0.009 when 
ignoring SN 2004dt). 

The overall trend and even the location of each object in 
the parameter space (given the measurement uncertainties) 
is quite similar to Figure [5j This is generally unsurprising 
given the results above involving the velocity gradient, since 
most (but not all) HV (normal-velocity) objects are also 
HVG (LVG) objects (Table|2|of [Wang et al. |2006l ISilverman 



et al. 2012d I. Furthermore, the Si II A6355 velocity at t = d 



(i.e., vq) can be approximated from the velocity gradient 
and, as expected, these velocities at maximum brightness 
correlate in almost exactly the same manner as what is seen 
in Figure [6] 

To investigate whether late-time observations of SNe la 
are related to the amount of circumstellar interaction a SN la 



undergoes, Forster et al. (20121 measured the EW of nar- 
row absorption from Na I D and Ca II H&K in the host 
galaxies of SNe la and compared them to nebular velocities. 
Using 212 BSNIP spectra of the 15 objects listed in Table[2] 
we attempt to measure the EW of host-galaxy Na I D and 
Ca II H&K absorption. Even in the relatively low-resolution 
spectra of the BSNIP sample, we are able to measure Na I D 
absorption in all 15 SNe la and Ca II H&K absorption in 
11 of them. The final EWs listed in the last two columns 
of Table [2] are the median of all EW measurements for each 
object (and the uncertainty is the standard deviation). 

There are 6 objects in the current work that were also 
studied by Forster et al.| ( |2012[ ) , and the EW measurements 
for all of these SNe la are consistent at the 2a level. In 
contrast to previous work ( Forster et al.|2012 |, the BSNIP 
data do not show a significant correlation between the EWs 



of Na I D and Ca II H&K absorption. Forster et al. (20121 



also found, at quite high significance, that objects with red- 
shifted nebular velocities have larger EWs of Na I D and 
Ca II H&K absorption. They interpret this result as indicat- 
ing that these objects have stronger circumstellar interaction 
or larger amounts of circumstellar material, and thus they 
claim that nebular velocities are directly related to SN la 
environments. 

The objects with redshifted ii no b in BSNIP tend to have 
stronger absorption from Na I D and Ca II H&K, but this 
is not always the case, and the difference between objects 
with blueshifted and redshifted nebular velocities is not sig- 
nificant. Therefore, we are unable to confirm the findings of 
|Forster et al.| ( [2012[ ) that the environments of SNe la directly 
affect nebular velocities (or vice versa). 

4-3.2 Photometric Comparisons 

As mentioned above, the photometric parameters used in 



this work can be found in Ganeshalingam et al. (2010). Of 



the 15 objects listed in Table [2] 14 have near-maximum- 
brightness light curves which allow for the measurement of 



Amis(B). We also analyse 17 additional objects from Maeda 
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Figure 7. Nebular velocity versus Amis (73). Shapes and colours 
of the data points are the same as in Figure[3] The six objects with 
Amis (7?) values which fall outside the "normal" range (0.9 < 
Ami5(/3) < 1.4 mag, denoted by the horizontal dashed lines) are 
labeled. The solid line is the linear fit to the data when ignoring 
these seven outliers. 



Figure 8. The deredshifted [Fe II] A7155 and [Ni II] A7378 
features of SNe 2006X, 1989M, and 1998bu, all of which have 
Ami5(73) RJ 1.1 mag. The solid vertical lines are the nebular ve- 
locity for both features of each object, while the dotted vertical 
lines are the rest wavelengths of the two features. 



et al. (2011 1 and Blondin et al. (20121. Most of them have 



0.9 < Amis(73) < 1.4 mag, which represents the range of 
"normal" SNe la (e.g., |Ganes halinga m et al.||2010| . This 
range is denoted by the horizontal dashed lines in Figure [7J 
where we plot w ne b against Amis (B) and have labeled the 
six objects that fall outside the normal range of Amis(73) 
values. 

When including all objects presented in the figure, the 
data are completely uncorrelated (Pearson correlation co- 
efficient of —0.009), similar to what has been seen pre- 
viously (Blondin et al. 2012 1. Even when the six outliers 



(SNe 1986G, 1999aa, 1999gh, 2003gs, 2003hv, and 2007on) 
are ignored, there is still no significant evidence of a correla- 
tion between ti ne b and Amis (B) (Pearson correlation coeffi- 
cient of —0.16). The solid line in Figure[7]is the linear fit to 
the data when ignoring these six outliers. Furthermore, the 
slope of this linear fit is consistent with zero, and it is clear 
that for a given Amis (B) there is a huge range of possi- 
ble v ne b values. This result, which matches that of previous 
studies (iMaeda et al.||2010a| iBlondin et al.||2012b, is illus- 



trated in Figure |8| where the [Fe II] A7155 and [Ni II] A7378 
features are plotted for three normal-luminosity objects all 
with Amis(B) « 1.1 mag. The solid vertical lines repre- 
sent the nebular velocity for both features of each object 
while the dotted vertical lines denote the rest wavelengths of 
the two features. Despite their nearly identical light curves, 
these three objects show significantly different nebular ve- 
locities with a range spanning ~4500 km s -1 . 
Silverman et al.| ( |2012c) ), (73 



In 



was defined as 



the observed B — V colour of the SN at 73-band maximum 
brightness. In Figure[9 we plot u nc b versus (73 — ^) max using 
this definition for the BSNIP dat43 as we ll as 14 additional 



SNe la from Maeda et al. (20111 and Blondin et al. (2012) 



In this parameter space there appears to be three outliers 



4 SNe 1989M and 1993Z do not have colour information available. 




-3-2-10 1 2 3 

Nebular Velocity (10 3 km s" 1 ) 

Figure 9. Nebular velocity versus observed (B — V) max . Shapes 
and colours of the data points are the same as in Figure [3] Two 
of the labeled outliers (SNe 1998bu and 2006X) have significant 
reddening from their host galaxies, while the third labeled outlier 
(SN 2003gs) shows little host-galaxy reddening. The solid line is 
the linear fit to all of the data. The circled points are objects in 
the low-extinction sample and the dotted line is the linear fit to 
this subset. 



(labeled in Figure|9|. Two of these, SNe 1998bu and 2006X, 
are known to have substantial extinction from their host 
galaxies (e.g., |Wang et al.||2008| |Matheson et al.||2008| ). On 
the other hand, the third outlier (SN 2003gs) shows no evi- 



dence for significant host-galaxy extinction ( Krisciunas et al. 
[2009] ). 

Utilising all data in Figure [9j « nc b is found to be cor- 
related (at the > 3<r level) with (B — V) (Pearson coef- 
ficient of 0.59), and the linear fit to these data is shown in 
the figure as the solid line. |Blondin et al. | ([2012 ) found no 
correlation between v ne b and intrinsic B — V colour at 73- 
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band maximum for their full sample or when removing the 
peculiar SN 2004dt. However, they see a correlation simi- 
lar to ours (with Pearson correlation coefficient of 0.82) for 
objects with v ae b > —2000 km s _1 ; note that most of the 
objects presented herein also have w M b > —2000 km s _1 . 

Since (B — V) max is an observed SN colour (and thus a 
mixture of intrinsic SN colour and reddening from the host 
galaxy), we also investigate possible correlations between 
Wnob and (B — V r ) max for various low-extinction subsets of 
our sample. Simply removing SN 2006X or both SNe 2006X 
and 1998bu from the analysis (since they are known to have 
large host-galaxy extinction) increases the Pearson coeffi- 
cient to 0.60 and 0.77, respectively (and both correlations 
are significant at the ~3<r level). Concentrating only on ob- 
jects with (B — V) max < 0.2 mag, we find a similarly strong 
correlation (Pearson coefficient of 0.64) again at the > 3<r 
level. 

On the other hand, a low-extinction sample was defined 
by Blondin et al. (20121, and no correlation between v ne b 
and intrinsic (B — V) max was detected. In this work we use 



similar criteria as Blondin et al. (20121 to define our low- 



extinction sample. A SN la is considered part of this sample 
if it resides in an E/S0 host galaxy, is found far (in projected 
distance) from a host galaxy of any type, or has a low value 
of A v according to MLCS2k2 ( |Jha et al.|2007[ ). Objects in the 
low-extinction sample are marked with a 'I' in Table |2| and 
are circled in Figure [9] In contrast to Blondin et al. (2012), 
we find a correlation between u nc b and observed (B — V) 
for the low-extinction sample (Pearson correlation coefficient 
of 0.63) at the ~2.5cr level, and the linear fit to this subset 
of SNe la is plotted in Figure [9] as a dotted line. 



5 CONCLUSIONS 

In this work we have analysed late-time (t > 100 d) opti- 
cal spectra of low-redshift (z < 0.1) SNe la which mostly 
come from the BSNIP dataset ( Silverman et al.|2012a |. The 
lone exception is SN 2011 by, and we present and analyse 
9 spectra of this object (2 of which are at late times). Af- 
ter removing objects that do not follow the Phillips relation 
(Phillips 1993) and spectra that do not cover either of the 
wavelength ranges of interest (namely the [Fe III] A4701, 
[Fe II] A7155, and [Ni II] A7378 features), we are left with 34 
SNe la with 60 spectra having t > 100 d. Of these, 7 were 



studied by Maeda et al. (2010a I, and those same 7 plus 2 



more were studied by Blondin et al. (20121. We add to the 



BSNIP data nebular spectral feature measurements of 20 



SNe la from previously published work ( Maeda et al.|[20TT 
Blondin et al.|20 12). This sample represents the largest set 
of late-time SN la spectra ever analysed. 

We attempted to measure the FWHM and velocity of 
the three nebular emission features mentioned above and 
are able to measure the [Fe III] A4701 line in all late-time 
data discussed in this work. On the other hand, we are only 
able to calculate velocities of the [Fe II] A7155 and [Ni II] 
A7378 features in 22 spectra of 15 SNe la. To successfully 
measure these features it seems that the SN la should be 
spectroscopically normal, the spectrum should have S/N > 
20 px _1 , and the spectrum should be older than ~160 d 



we calculate v ne h from the mean of the individual velocities 
of the [Fe II] A7155 and [Ni II] A7378 features. 

The FWHM of [Fe III] A4701 is anticorrelated with 
Am 1 c,(B). However, if the 6 objects with Amis(B) > 
1.6 mag are removed, the correlation effectively disappears. 
The same results are found when using only spectra from 
250-300 d past maximum brightness. Thus, this possible 
relationship seems to be driven solely by the most under- 
luminous SNe la. This feature is blueshifted by ~ 2000- 
4000 km s _1 at t = 100 d, and the amount of blueshift 
decreases approximately linearly with time. However, this 
temporal evolution is quite slow, with a typical change in 
velocity of only ~10-20 km s _1 d _1 . 

The nebular velocity likewise does not change much 
with time (typical values are a few km s _1 d _1 ). Further- 
more, the BSNIP data weakly support the finding ofjMaeda] 
et al. (2010a I that « ne b and v are correlated. Nearly all HVG 



objects have redshifted nebular lines and most LVG objects 
have blueshifted nebular spectral features. A similar result 
is found when using the near-maximum-brightness velocity 
instead of the velocity gradient and assuming HV (normal- 
velocity) objects are also HVG (LVG) objects (which is a 
reasonable assumption, e.g., Silverman et al. 2012d I . The 



data studied in this work also show no significant correla- 
tion between the EWs of Na I D or Ca II H&K absorption 
and Wnob! thus, we are unable to confirm previous claims that 
the environments of SNe la directly affect nebular velocities, 
or vice versa ( Forster et al.|20T2 k No significant correlation 
is found between v ne b and Amis(-B), and SNe la with nearly 
identical light curves can have a wide range of nebular veloc- 
ities, spanning ~4500 km s _1 in Figure|8l Our data also in- 
dicate a correlation between observed (B — V) an d iWb- 
Although the combined dataset analysed in this work 
constitutes the largest late-time SN la spectral dataset ever 
studied, it still contains only a handful of objects. The anal- 
ysis herein could certainly be improved and extended by 
obtaining more nebular SN la spectra. While BSNIP (and 
the CfA, Blondin et al.| |2012) has gathered, published, and 
analysed their SN la spectral data, larger transient searches 



such as the Palomar Transient Factory (PTF; Rau et al. 
2009"}|Lawet al.|2009| > and Pan-STARRS ( |Kaiser et al.|2002 \ 



past maximum brightness. Following Maeda et al. (2010a I, 



will be able to add significantly to the number of late-time 
SN la spectra. 
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Figure Al. BVRI light curves of SN 2011by from KAIT. 
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APPENDIX A: OBSERVATIONS OF SN 2011BY 



SN 2011by was discovered by |Jin fc Gao| ( |2011[ ) on 
2011 Apr. 26.8 UT in NGC 3972 with J2000.0 coordinates 
a = ll h 55 m 45.56 s , 5 = +55°19'33'.'8, and it was classified 



-4 d after discovery as a young SN la by Zhang et al. ( 2011 1 



Immler & Russell (2011 I detected X-ray emission coincident 
with SN 2011by ~40 d after discovery using Swift/XKT, but 
the emission was found to be unassociated with the SN using 
Chandra observations taken ~3 weeks later ( Pooley|2011[ ). 

Broad-band BVRI photometry of SN 2011by was ob- 
tained using the 0.76-m Katzman Automatic Imaging Tele- 



scope (KAIT) at Lick Observatory (Li et al. 2000 Filip 



penko et al. 



The data 



2001 1 and is shown in Figure Al 
were reduced using standard procedures, the details of which 

poTol. The optical 



light curves of SN 2011by indicate that it reached B-band 
maximum brightness on 2011 May 9.9 UT and that it had 
Amis(B) = 1.14 ± 0.03 mag. 

We obtained 7 near-maximum-brightness spectra of 
SN 2011by using the Lick/Kast double spectrograph. They 
span ~12 d before to ~15 d after maximum brightness, as 
shown in the top panel of Figure |A2| In addition, we used 
Keck/LRIS to acquire 2 late-time spectra at 206 and 310 d 
past maximum (shown in the bottom panel of Figure |A2[ ) , 
the spectral feature measurements of which can be found at 



Figure A2. Near-maximum-brightness (top) and late-time (bot- 
tom) spectra of SN 2011by, labeled with age relative to the time 
of maximum brightness. The data have had their host-galaxy re- 
cession velocity removed and have been corrected for Galactic 
reddening. 



the end of Table[T] Details of all of our spectral observations 
of SN 201 lby can be found in Table |AT) 

Consistent with the light curves, the spectra indicate 
that this object is quite normal; the SuperNova IDentifica- 



tion code (SNID; Blondin & Tonry 20071, as implemented 
by [Silverman et al.| ( |2012a[ ), was run on all of our spectra 
of SN 2011by and it is definitively classified as spectroscop- 



ically normal. As in Silverman et al. (2012dl, we calculated 



the near-maximum-brightness velocity based on the mini- 
mum of the Si II A6355 absorption feature and find that 
it shows a "normal" velocity near maximum brightness, as 
opposed to some SNe la which exhibit higher-than-average 
expansion velocities (e.g., Wang et al. 20091. Furthermore, 



since we have multiple spectra near maximum brightness, 
we can calculate the velocity gradient (e.g., |Benetti et al.| 
20051); SN 201 lby falls squarely in the LVG group. 



Individual objects (some found by amateur astronomers 
and some by large-scale surveys) discovered in nearby galax- 
ies soon after explosion will be invaluable resources for ex- 
tending our understanding of SNe la. SN 201 lby represents 
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UT Date 


Age a 


Range (A) 


Airmass b 


Exp (s) 


2011 Apr. 28.2 


-11.6 


3440- 


-10200 


1.05 


600 


2011 May 4.2 


-5.6 


3440- 


-10200 


1.06 


600 


2011 May 6.2 


-3.7 


3440- 


-10200 


1.05 


600 


2011 May 10.2 


0.3 


3438- 


40196 


1.06 


600 


2011 May 14.2 


4.3 


3456- 


-10200 


1.07 


600 


2011 May 22.3 


12.3 


3464- 


-10200 


1.15 


600 


2011 May 25.3 


15.4 


3436- 


-10200 


1.27 


600 


2011 Dec. 2.6 C 


206 


3400- 


-10138 


1.50 


450 


2012 Mar. 15.5 C 


310 


3400- 


-10200 


1.32 


900 



a Rest-frame days relative to the date of £?-band maximum 
brightness, 2011 May 9.9 (see Section 
b Airmass at midpoint of exposure. 



2.1 1 



c These observations used LRIS ( |Oke et al.|1995| on the 10-m 
Keck I telescope. The others used the Kast spectrograph on 
the Lick 3-m Shane telescope I Miller fe Stone|1993 i. 



an exquisite case study for the current work. In recent years 
there have been a handful of other nearby SNe la found 



soon after explosion (e.g., Nugent et al. 2011 Foley et al 



2012 Silverman et al. 2012b) that will also serve as excel- 



lent individual case studies once late-time spectra have been 
obtained and analysed. 
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